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Abstract 

The static dipole polarizability of noble gas atoms confined by impenetrable spheres and spherical 
shells is studied using the B-spline random phase with exchange approximation. The general trend 
in dipole polarizabilities across the noble gas sequence shows a decrease in the dipole polarizability 
as the volume of the confining impenetrable sphere is reduced and a large increase in the dipole 
polarizability for confinement by impenetrable spherical shells as the inner shell radius is increased. 
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I. INTRODUCTION 


The study of atoms under confinement is a subject of long standing interest, for a selection 
of reviews see |l]- jo] . For many unconfined atoms, their polarizabilities are known to a 
high degree of accuracy, see {?J and js]. An interesting question arises as to the effect of 
confinement on the polarizability of an atom. 

The majority of studies investigating the polarizabilities of atoms under confinement have 
been carried out for hydrogen, see [9|-[ll| for some early work and 12]- 18] for more recent 


19] has been followed by density functional 


calculations. For helium, an early calculation 

20], configuration interaction ]2l|, coupled cluster |22] and variational 


23 


confining potentials studied including cylindrical ]24] and spheroidal [25|. For heavier sys- 


calculations with 


terns, one-electron model potentials have been used to study the static dipo 


of alkali metals under hard wall confinement 


26] and debye screening 27]- 


e polarizability 


29]. The static 


dipole polarizability and hyperpolarizability of beryllium was studied under harmonic, rect¬ 
angular and gaussian confining potentials using a basis of explicitly correlated gaussians jsoj]. 
For the polarizabilities of noble gas atoms under hard wall confinement, an early calculation 
studied the confined argon atom in a Thomas-Fermi mode 


32]- 34] and a time-dependent density-functional theory 


311] . Hartree-Fock calculations 


35] calculation have since been 


performed along the noble gas sequence. The polarizabilities of confined noble gases are of 
interest in interpreting the density dependence of atomic polarizabilities 36] and studies of 
noble gas endofullerenes 37]-[39]. 

Here, non-relativistic random phase approximation with exchange (RPAE) calculations 
of the static dipole polarizability are performed for noble gas atoms under hard-wall confine¬ 
ment. The RPAE method 0-1 I 42 I has been shown to give dipole polarizabilities accurate to 
within a few percent for the noble gases when compared to more sophisticated calculations 


43], with relativistic effects having a small effect on polarizabilities for the noble gases 44], 


The rest of the paper is structured as follows: Section II details the RPAE approximation 
to the calculation of the static dipole polarizability and its numerical implemention using 
B-splines. In section III A, results are presented for the static dipole polarizabilities of 
hydrogen and noble gas atoms confined by an impenetrable sphere. Section III B investigates 
the case of a spherical shell potential. Finally, Section IV concludes with a summary. Unless 
otherwise stated, atomic units are used throughout. 
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II. THEORY 


A. Calculation of the dipole polarizability 


For a hydrogen atom in the ground state, the static dipole polarizability can be calculated 
via the relation, 


a 


= 2E 


|(np|z|ls)p 


( 1 ) 


n ‘-rip £ls 

For noble gas atoms, electron correlation beyond the Hartree-Fock approximation plays 
an important role and must be accounted for. The dipole polarizability for noble gas atoms 
will be calculated via the RPAE approximation 0-i S i . In the following equations for the 
RPAE vertex matrix element, the summation is over hole states, or states below the Fermi 
level (< F), and excited electron states above the Fermi level ( > F). 

For (u > F) and (/j, < F), 

( \ 

E-E 


(u\A(u)\n) = (v\d\fi) 
and for (u < F), (// > F), 
(v\A(u)\n) = (v\d\y) 


u'>F 




y<F v'< 


fJ 


00 — Ei/i T Ep,' + z5(l — 2Uj/) 


/ 




E-E 

v'>F [i'>F 

V'<F v'<Fj 


(W\vWn)-W*\vWn)), AAy)W ' (3) 


00 — £ v > + Epi + f<5(l — 2 n u >) 
where id is an infinitesimal imaginary positive, and 


n,/ = 


0 for v' > F 
1 for v' < F 


( 4 ) 


Numerically, equations (J2J) and d2J) are a set of linear algebraic equations, which can be 
written in block matrix form, 
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\y 
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d 
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/ 


fd r 2 
R 2 Vi 


Xi 0 
0 X2 
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X 


\y 


( 5 ) 


where Xi and y 2 are energy denominators, x and y represent (v\A{oo)\y) for v > F, y < F 
and v < F , /x > F respectively, d are the Hartree-Fock dipole matrix elements and Vi, V 2 
are Coulomb matrices. 
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Once this equation has been solved for the RPA vertex, the dipole polarizability can be 
calculated via the relation: 


/ n = y (fJ>\d\i')(i'\A{u)\fj) + y (n\d\u)(u\A(u)\fi) 
“7 U - e v + e^ + id Z^ F -u + e„ - + i6 

fi<F /i>F 

In this work the static dipole polarizability, uj = 0, is calculated. 


( 6 ) 


B. Numerical Implementation 


Vsphiy ) 


(7) 


V she u(r ) = < 


The Hartree-Fock equations will be solved for noble gas atoms confined by a potential. 
Two forms of the confining potential will be compared here. First, an impenetrable spherical 
potential of the form, 

0 for r < r c 
oo for r > r c 

where r c is the radius of the confining potential. Secondly, an impenetrable spherical shell 
potential of the form. 

oo for 0 < r < r i n 

0 for r in < r < r ou t ( 8 ) 

_ oo for r > r 0 ut 

where r in is the inner radius and r out is the outer radius of the spherical shell potential. As 
the radius of the confining potential may lie within the radius of the unconfined atom, it is 
necessary to incorporate potentials ((71) and ((8]) into the numerical solution of the Hartree- 
Fock equations 45|, 46]. For all noble gas atoms, the Hartree-Fock equations were solved 
using a cut off radius of 30 a.u. and 6000 integration points. Note that for heavier atoms, 
convergence of the Hartree-Fock equations for small radii of the confining potential is more 
problematic. Results will be presented in this work down to the minimum radius of the 
confining potential for which convergence was obtained. 

As polarizabilities of atoms can involve substantial contributions from continuum states, 
an efficient means of spanning the single particle energy continuum is needed. Here, B-spline 
basis sets will be used [f7, 48]. B-splines are piecewise polynomials defined on a particular 
knot sequence. For confinement by an impenetrable spherical potential, an exponential knot 
sequence is used of the form, 

h = t 2 = ■ ■ ■ = 4 = 0, (9) 
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( 10 ) 


tn+\ tn+2 ' tn+k f c . 


For t k+1 ^t n , 


U = r 0 (exp (P(i - k)) - 1) 

M* + 1) 


P = 


( 11 ) 

( 12 ) 


(n + 1 — k) ’ 

where is a parameter which affects the smallest spacing of the knot sequence. For con¬ 
finement by an impenetrable shell potential, an exponential knot sequence is used of the 
form, 


For t k +i->t n , 


t\ — ^2 — ' ' 


(13) 

^n+1 ^n+2 

t'n+k ^out* 

(14) 

U = r in + r o (exp(/3((f - k)) - 1 ) , 

(15) 

l n ( (rout-rm) + ^ 

a K r 0 

P (n + 1 - k) ' 

(16) 


A value of = 0.001, with n — 60 splines of order k = 9 was used in this work. 
By expanding the radial wavefunctions in terms of B-splines, 

flW = E4' ) aw, 


(17) 


and substituting into the atomic hamiltonian for hydrogen or the Hartree-Fock equation for 
noble gas atoms, 

H^P l {r)=eP l {r\ (18) 

a generalised eigenvalue problem is obtained for angular momentum l that can be solved to 
obtain a complete set of states suitable for the calculation of the atomic polarizabilities via 
equations (|T]) - (|6lh The present calculations were carried out up to a maximum orbital angular 
momentum of l = 10. For confinement by an impenetrable spherical sphere, the boundary 
conditions P/(0) = Pi{r c ) = 0 are implemented by removing the first and last splines, leaving 
n-2 electron eigenstates for each orbital angular momentum l. Similarly, for confinement by 
an impenetrable spherical shell, the boundary conditions P/(r in ) = P/(r out ) = 0 are also 
implemented by removing the first and last splines, leaving n-2 electron eigenstates for each 
orbital angular momentum l. For more details on the numerical implementation of the 


RPAE approximation, see chapter 2 of 


42 ] 
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III. RESULTS 


A. Static dipole polarizabilities of atoms confined by an impenetrable spherical 
potential 


In order to benchmark the B-splinc basis used in this work, equation © is used to 
calculate the static dipole polarizability of a hydrogen atom confined in an impenetrable 
sphere, with the results tabulated in table [TJ Comparison with the results of perturbation- 


Hyllcraas 


18] and perturbation-numerical 


14] calculations shows excellent agreement with 


the B-spline results. Comparison with the lower bound Kirkwood 18] and Buckingham 18] 


results shows that the higher order terms included in the Buckingham approximation 


50] 


eads to a more accurate lower bound on the polarizability than the Kirkwood approximation 


5111 - 


table I: Static dipole polarizabilities for a hydrogen atom confined in an impenetrable sphere. 
R Present results Kirkwood [18] Buckingham [18] Hylleraas [18] Perturbation m 


2.0 0.34255811 

0.3401420 

0.3401565 

0.3425581 

0.342558 

4.0 2.37798233 

2.2908060 

2.3578222 

2.3779823 

2.377982 

6.0 4.05814049 

3.7054757 

4.0471005 

4.0581405 

4.058140 

8.0 4.45396473 

3.9750587 

4.4527846 

4.4539647 

4.453965 


To calculate the polarizability for multi-electron atoms via the RPAE approximation, 
bound states found from solving the Hartree-Fock equations are used to construct a complete 
set of states in a B-splinc basis. Tables UTl and II I II compares the present Hatree-Fock subshell 
energies for helium and neon with the Roothaan-Hartree-Fock results of { 52 1. For the energy 
of the Is subshell of helium there is good agreement between the calculations, while for the 
Is, 2s and 2p subshells of neon, the present subshcll energies are more deeply bound than 
those of 52], To further investigate this discrepancy, Roothaan-Hartree-Fock calculations 
for Neon similar to those of 52] could be performed, examining convergence of the subshell 
energies with respect to both the orbital basis set size and the optimisation of the basis set 
exponents. 

Table [TV] contains static dipole polarizabilities calculated in the RPAE approximation for 
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TABLE II: Hartree-Fock energies for the Is subshell of He confined in an impenetrable sphere. 


R Eis (present work) Ei s [52] 
2.0 -0.65754880 -0.6593 

4.0 -0.91252115 -0.9125 

6.0 -0.91787985 -0.9174 


TABLE III: Hartree-fock energies for the Is, 2s and 2p subshells of Ne confined in an impenetrable 
sphere. 


R Is (present work) Ei s [52] E 2 S (present work) E 2 S [52] E 2 P (present work) E 2 P [52] 


2.5 -32.41538078 

-32.3831 -1.72747430 

-1.7084 -0.62658124 

-0.6076 

3.0 -32.63113594 

-32.5892 -1.84794644 

-1.8239 -0.76329739 

-0.7386 

3.5 -32.71673490 

-32.6763 -1.89664949 

-1.8731 -0.81625604 

-0.7924 

4.0 -32.75109004 

-32.7166 -1.91699524 

-1.8963 -0.83731963 

-0.8169 


noble gas atoms confined in an impenetrable sphere in comparison to unconfined polariz¬ 
abilities, with the results displayed graphically in figure CD It is seen that as the confining 


potential radius increases, the do 


arizability converges smoothly to closely agree with the 
441] , demonstrating the small influence on the polarizabilities 


relativistic RPA calculations of 
of relativistic effects. The non-relativistic RPA calculations of [40| are slightly lower than 


the present results with the exception of xenon where the results of [40] are hig 


rer. The 


differences are within the estimated 5% uncertainty of the earlier calculations 40]. As the 


confining potential radius decreases, the polarizability decreases monotonically. Compari¬ 
son with more precise hylleraas calculations for helium 53] and relativistic coupled cluster 
calculations for neon, argon, krypton and xenon 43] shows that the non-relativisitic RPAE 
method is able to recover the static dipole polarizability to within a few percent. 

In figure [2l the present RPAE polarizabilities for helium are compared to Hartree-Fock 
(HF) 32], density functional theory (DFT) 20] and time-dependent density-functional the¬ 
ory (TDDFT) 35] calculations. The Hartree-Fock polarizabilities of jssj] are larger than 
the current results and the other calculations at all box radii. The time-dependent density- 
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TABLE IV: Static dipole polarizabilities calculated in the RPAE approximation for nobles gas 
atoms confined in an impenetrable sphere. 


R 

He Ne 

Ar 

Kr 

Xe 

1.8 

0.27149963 




2.0 

0.35985693 




2.5 

0.60574486 0.99099895 




3.0 

0.84391160 1.39793917 

3.96443519 



3.5 

1.03453905 1.74913386 

5.36273398 

6.64086248 


4.0 

1.16491541 2.00968621 

6.72840349 

8.67526518 

11.34371382 

5.0 

1.28514057 2.27781370 

8.89416026 

12.30276614 

17.33513804 

6.0 

1.31534685 2.35591238 

10.07330554 

14.67047514 

22.05689052 

7.0 

1.32113843 2.37303121 

10.54979313 

15.82775188 

24.91468796 

8.0 

1.32206986 2.37613351 

10.70305446 

16.27562844 

26.28682580 

9.0 

1.32220435 2.37663254 

10.74468581 

16.42077761 

26.83277612 

10.0 

1.32222378 2.37671053 

10.75476508 

16.46192981 

27.02021394 

15.0 

1.32223132 2.37673811 

10.75782869 

16.47631699 

27.10088802 

20.0 

1.32223265 2.37674263 

10.75792504 

16.47654554 

27.10156072 

30.0 

1.32223336 2.37674498 

10.75797675 

16.47666561 

27.10188992 

RPAE [40] 

2.30 

10.73 

16.18 

27.98 

RRPA[44] 

1.322 2.38 

10.77 

16.47 

26.97 

RCCSDT[43] 

2.697 

11.22 

16.80 

27.06 

Hylleraas [53] 1.38376079 





functional theory polarizabilities of 35] are higher than the current results above a confining 
potential radius of 3 a.u. and fall off more strongly below a confining potential radius of 
3 a.u. Excellent agreement with the present results is seen with the exact exchange, EXX, 
density functional results of 20]. The exchange-only local-density approximation, XO-LDA, 
and exchange-correlation local-density approximation, XC-LDA, polarizabilities trend higher 
than the present results for larger confining potential radii, but agree better for smaller radii. 

For Neon, figure [3] compares the present RPAE polarizabilities to Hartree-Fock 32] and 
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FIG. 1: Static dipole polarizabilities of the noble gas sequence as a function of the confining 
potential radius. Solid line with circles, helium; dotted line with squares, neon; dashed line with 
diamonds, argon; dot-dashed line with triangles, krypton; dot-dashed-dashed line with pluses, 
xenon. 


time-dependent density-functional theory 35] calculations. Similarly to the case of helium, 


the Hartree-Fock polarizabilities of 


32] are substantially higher than the current results. 


The time-dependent density-functional theory polarizabiltiies of |35j are in good agreement 
with the present results, with a distinctive sigmoidal shape seen in the polarizabilities as a 
function of the confining potential radius. 
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FIG. 2: Static dipole polarizabilities of helium as a function of the confining potential radius. Solid 


line with circles, present results; dotted line, HF 


TDDFT 


351 ] (digitized using g3data 


line, DFT (XO-LDA) 


MD; 


321 ] (digitized using g3data 


dot-dashed line, DFT (EXX) 


20|; dot-dot-dashed, DFT (XC-LDA) 


541]): dashed line, 


20(; dot-dashed-dashed 


2C|. 
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FIG. 3: Static dipole polarizabilities of neon as a function of the confining potential radius. Solid 
line with circles, present results; dotted line, 


(digitized using g3data 


321 ] (digitized using g3data 


54]]); dashed line, 


35] 


m 


ii 











B. Static dipole polarizabilities of atoms confined by an impenetrable spherical 


shell potential 


It has been demonstrated that confinement of noble gas atoms by an impenetrable spher¬ 
ical potential leads to a monotonic decrease in the polarizability as the radius, and hence 
volume of the confining potential is reduced. A natural extension is to an impenetrable 
spherical shell potential, c.f. equation (JHJ) - This potential may be of interest to the case 


of an atom confined by a negatively charged fullerene shell 


551 ]. The polariza 


jility of hy¬ 


drogen confined in an impenetrable shell potential has been studied in 56]- 60]. Here, the 
confinement of helium and neon in a shell potential will be examined. 

As for the case of an impenetrable spherical potential, the present calculations are first 
benchmarked for hydorgen. Table |V] compares the static dipole polarizabilities calculated 
using equation ([I]) against existing calculations 56] for hydrogen confined in a spherical 
shell potential of outer radius 10 a.u. with a varying inner shell radius. There is seen to 
be close accord between the present results and those obtained in [56] using a pertubation 
numerical method and the Buckingham approximation 50] • Similarly to confinement by an 
impenetrable spherical potential, results obtained using the Kirkwoord approximation 51] 
tend to underestimate the polarizabilities. The polarizabilities increase rapidly as the inner 
shell radius is increased. 


TABLE V: Static dipole polarizabilities for hydrogen confined by an impenetrable spherical shell. 


r-m r ou t Present results Perturbation [56] Kirkwood [56] Buckingham [56] 

2.0 10.0 485.28323894 474.3865 420.8943 465.6254 

4.0 10.0 1393.53470841 1372.0345 1043.0354 1352.7657 


RPAE calculations of the polarizabilities of helium and neon confined by a spherical shell 
potential of outer radius 10 a.u. with a varying inner shell radius are presented in table IVTl 
The trend in the polarizabilities for helium and neon is similar to the case of hydrogen, with 
the polarizability for neon rising more strongly as the inner shell radius is increased than 
that of helium. 

The increase of the polarizability of neon as the inner shell radius is increased can be 
explored by examining the subshell energies, see table IVIII It is seen that a subshcll re- 
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TABLE VI: Static dipole polarizabilities calculated in the RPAE approximation for helium and 
neon confined by an impenetrable spherical shell. 


Tin r ou t He Ne 


0.0 

0.25 

0.5 

0.75 

1.0 

1.25 

1.5 
1.75 
2.0 

2.25 

2.5 


10.0 1.32222378 
10.0 12.16238797 
10.0 31.94196072 
10.0 60.45733816 
10.0 96.27102702 
10.0 137.37472602 
10.0 181.62662627 
10.0 227.70588636 
10.0 274.39609861 
10.0 320.86226234 
10.0 366.59086983 


2.37671053 

73.25472372 

164.42341795 

251.97880982 

335.74822122 


arrangement 61] occurs as the inner shell radius is increased with the 2s subshell becoming 
more weakly bound than than the 2p subshcll as electrons are excluded from penetrating 
into the region near to the nucleus. As the inner shell radius is increased the binding of the 
2s orbital is decreased leading to a large increase in the polarizability. 


TABLE VII: Hartree-fock energies for the Is, 2s and 2p subshells of neon confined by an impene¬ 
trable spherical shell. 


Hri r ou t Els 


E2s 


E 


2p 


0.0 10.0 -32.77243555 -1.93039045 -0.85040971 
0.05 10.0 -11.67496109 -0.59126878 -1.30650997 
0.25 10.0 -3.72928450 -0.24957000 -1.52410597 
0.5 10.0 -2.18429623 -0.18295290 -1.23595070 
0.75 10.0 -1.56730085 -0.14556270 -1.00631972 
1.0 10.0 -1.22548285 -0.11740385 -0.84030892 
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IV. CONCLUSIONS 


The static dipole polarizabilities of noble gas atoms confined in impenetrable spherical 
spheres and impenetrable spherical shell potentials has been studied using the RPAE ap¬ 
proximation. Confinement in an impenetrable sphere leads to a decrease in the polarizability 
as the radius of the sphere is decreased. However confinement by an impenetrable spher¬ 
ical shell potential leads to a large increase in the polarizability as the inner shell radius 
is increased. This offers the potential of tuning the polarizability of an atomic system by 
confinement. Topics in the polarizabilities of confined multi-electron atoms that require fur¬ 
ther investigation include the influence of relativistic effects 62] and exploration of a range 
of confining potentials beyond the hard-wall confinement approximation, such as repulsive 


63j, 64] or attractive {t?]] penetrable spherical shell potentials, with a view to exploring both 


higher multipole polarizabilities and dynamic polarizabilities. 
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